Small-scale properties of the KPZ equation and dynamical symmetry breaking 
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A functional integral technique is used to study the ultravi- 
olet or short distance properties of the Kardar-Parisi-Zhang 
(KPZ) equation with white Gaussian noise. We apply this 
technique to calculate the one-loop effective potential for the 
KPZ equation. The effective potential is (at least) one-loop 
ultraviolet renormalizable in 1, 2, and 3 space dimensions, 
but non-renormalizable in 4 or higher space dimensions. This 
potential is intimately related to the probability distribution 
function (PDF) for the spacetime averaged field. For the re- 
stricted class of field configurations considered here, the KPZ 
equation exhibits dynamical symmetry breaking (DSB) via 
an analog of the Coleman- Weinberg mechanism in 1 and 2 
space dimensions, but not in 3 space dimensions. 
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In understanding the onset of spatio-temporal pattern 
formation in systems out of equilibrium, it has proven 
extremely useful to begin the pattern formation and se- 
lection analysis by first obtaining and then classifying 
all the static and spatially homogeneous states allowed 
by the time dependent partial differential equations that 
model the system in question [0. In this way, one can 
decide whether the system will exhibit Hopf bifurcations 
and/or Turing instabilities and get a handle on the qual- 
itative nature of the patterns expected to emerge. This 
can be followed up by an amplitude analysis of the fluc- 
tuations about these static and homogeneous configura- 
tions. The unstable modes are the ones that lead to non- 
trivial patterns. For out-of-cquilibrium systems coupled 
to noisy environments (or with inherent internal noise) it 
is important to know how the stochastic sources can alter 
and shift these static and homogeneous configurations, 
since these affect the onset of the pattern-forming insta- 
bilities. Here we show how the effects of noise on these 
configurations can be computed using functional integral 
methods, and reveal an intrinsically ultraviolet one-loop 
phenomenon not captured by the more common meth- 
ods of analysis described above: namely the occurrence 
of dimension dependent dynamical symmetry breaking. 
In this paper we focus attention on the Kardar-Parisi- 
Zhang (KPZ) equation 



(1) 



The KPZ equation arises in a host of seemingly dis- 
tinct physical contexts and problems ranging from mod- 
els of turbulence |^,p|, interface growth, driven diffu- 
sion and flame fronts" |] , directed polymers in a random 
medium , certain lattice gases with hard-core exclu- 
sion [H, and even structure development in the early uni- 
verse^. In the fluid dynamics interpretation (when the 
KPZ equation is used as a model of Burgers turbulence), 
the fluid velocity is u = —Vcf) and the KPZ field plays the 
role of a velocity potential. In the surface growth inter- 
pretation 0(2;, t) is the local height of the surface, defined 
over a two-dimensional plane |^ . The constant "tadpole" 
term i^o in (|l|) is necessary for the ultraviolet regulariza- 
tion of the KPZ equation^ After renormalization, 
we will argue that the tadpole can and should be set to 
zero. 

If the noise 77 is Gaussian, all stochastic averages are 
encoded in the generating functional^ 



Z[J]^ / (Vcj)) exp 



exp 



J0 



(2) 



In general, a functional Jacobian determinant must be 
included. For the KPZ equation this determinant is a 
field-independent constant [lO| . For translation- invariant 
Gaussian noise, we split its two-point function into an 
amplitude, A, and a shape, g2{x,y) 



, dcf 



with the convention that 



(3) 



^In this regard it is interesting to point out that small scale 
properties of a randomly stirred fluid were studied some years 
ago by Yakhot, who found that the ultraviolet renormalizabil- 
ity of the forced Navier-Stokes (NS) equation required the 
addition of a term not originally present in the bare NS equa- 
tion 0. 

^This "direct" functional integral formalism is different (yet 
equivalent) from the more extended MSR approach p^ . 
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J g^\x, t) d^x dt = I = ~g^\k = 0, = 0). (4) 

The amplitude A is the loop-counting parameter for this 
theory @. 

There are two important symmetries of the KPZ equa- 
tion that are relevant for our analysis. First, we have 



0^0 + c(t), 
dc{t) 



Fq Fq 



dt 



(5) 
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In the fluid dynamics interpretation this symmetry is a 
"gauge transformation" of the scalar field that does not 
change the fluid velocity v. In the surface growth inter- 
pretation this symmetry is a (Type I) Galilean transfor- 
mation. It can be used to eliminate any spurious mo- 
tion of the background field. The second symmetry is 

(k«i|) 



X ' ^ X — X — A 6 

(f , t) (l)'{x', t') = , t) - e • X. 



(7) 
(8) 
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In the fluid dynamics interpretation this symmetry is 
equivalent to a (Type II) Galilean transformation of the 
fluid velocities 

v^v'::^v + e. (10) 

In the surface growth interpretation this symmetry 
amounts to tilting the coordinate system at an angle to 
the vertical, with 



tan(0)-||e|| 



(11) 



This Type II Galilean transformation is an exact invari- 
ance of the zero-noise KPZ equation, but once noise is 
added to the system, it will remain a symmetry only if 
the noise is translation-invariant and temporally white. 

A few words are in order regarding the choice of field 
configurations used for calculating the effective poten- 
tial. In quantum field theory (QFT) patterns of spon- 
taneous and dynamical symmetry breaking are revealed 
using static and homogeneous fields, since these configu- 
rations minimize the positive kinetic energy contribution 
to the total energy. The QFT ground states (basic field 
configurations) are therefore determined by searching for 
the minima of the effective potential. It turns out that an 
effective potential for the KPZ equation can be defined 
for a static, but spatially inhomogeneous field configura- 
tion, if it has the form 

(/, = -w-f. (12) 

In the hydrodynamic interpretation this corresponds to 
a stationary and homogeneous fluid flow Q with velocity 
V. In the surface growth interpretation ||w|| corresponds 
to a constant slope of the surface |0] (linear ramp). Due 



to the structure of the stochastic field theory (g) lead- 
ing to the KPZ equation, this is the most general static 
configuration allowing one to define an effective poten- 
tial. This is because the effective potential must be a 
function of a static field that satisfies the heat equation: 
dt4>— v'S/'^4> = 0. Any other field configuration more gen- 
eral than ( p^ ) (and/or time dependent) would force us to 
deal with the complexities of the effective action. 

Following the general construction in |jl^, the KPZ 
zero-loop effective potential is given by 



v.. 



zero — loop [ 



{Fo + iXv^y 



F' 
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The effective potential is generally calculated in terms of 
a reference background field (f>Q = —vq ■ x. We have used 
the Type II symmetry to set vq = 0. This potential is 
formally equivalent to that of At/)^ QFT - with the veloc- 
ity V playing the role of the quantum field 0qft- Even at 
zero-loops we see that for Fq < 0, the effective potential 
takes on the "Mexican hat" form, so that the onset of 
spontaneous symmetry breaking (SSB) would not be at 
all unexpected |lj,|l^. However, because of the Type I 
symmetry, the renormalized value of Fq must be set to 
zero, and SSB is not encountered in the KPZ equation. 
(There is an instructive analogy with the renormalization 
program for massless A(/)^ QFT. It is well known that for 
this QFT one still needs to introduce a bare mass param- 
eter that must be kept for all intermediate stages of the 
regularization. It is only after renormalization is com- 
plete that it makes sense to set the renormalized value 
of the mass to zero.) Considerably more subtle is the 
onset of DSB which can be detected only via a one-loop 
computation. 

Evaluating the general expression for the one-loop ef- 
fective potential computed in for the case of the KPZ 
equation, we obtain 



i(Fo + iA«2)2 



{2n) 
O) + 0{A^) 



~g2ik,iu)X{Fo + iXv^)e 



{uj-Xv-kf 



l^2(/c2)2 

(14) 



Because we have chosen noise which is white in time, 
so that g2{k,uj) — > .92(^)7 we can shift the integration 
variable from oj to uj — Xv ■ k. The resulting integral 
becomes 

VH = i [{Fo + iXvr - FS] +iAl 

'uj^ + i^^iPf + g2(fc)A(Fo + ^Xv^)P 



X loa 



LO^ + iy^{k^)'^ +g2{k)XFok'^ 



0(^2 )_ 



We point out that the denominator of the logarithm in 
the integrand is independent of the velocity field v. This 
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important result is a consequence of the fact that the 
Jacobian functional determinant encountered in is 
constant for the KPZ equation. One can perform the 
frequency integral exactly to obtain 



(15) 



The bare KPZ potential ( p^ ) contains terms proportional 



to , , and v . (The one-loop contribution (|15|), ex- 
panded in powers of , has terms proportional to v"^" .) 
Let us now take the spatial noise spectrum to be cutoff 
white noise, i.e., 



ff2(fc) = 52(||fc||) =e(A-A). 



(16) 



With this choice of noise, the term is proportional to 
A*^, the v'^ term to h.'^~'^ , and the term to A''^''. In or- 
der to be able to absorb the infinities into the bare action 
we must have d < A. Note that d = 4 is a marginal case, 
and one might suspect that d = 4 is a critical dimension 
for the KPZ equation That is: the KPZ equation 

(subject to white noise) is one-loop ultraviolet renormal- 
izable only in 1, 2, and 3 spatial dimensions. Even so, 
one-loop renormalizability requires an explicit tadpole. 
(Without a tadpole there is no term proportional to 
in the zero-loop potential, and hence no possibility of reg- 
ularizing the leading divergence.) Strictly speaking, the 
claim of one-loop renormalizability requires investigation 
of the effective action (wave- function renormalization). 
A one-loop ultraviolet renormalization of the effective ac- 
tion, valid for time dependent and inhomogeneous fields, 
has been carried out which fully supports the effective 
potential calculations presented here |1^ ]. 

The ultraviolet renormalizability of the KPZ equation 
depends critically on the ultraviolet behavior of the noise. 
Let us suppose that the noise is power-law distributed in 
the ultraviolet region with g2{k) w {ko/ky 9(A — k). 
In this case the nth term in the expansion has ultravio- 
let behavior proportional to y-^r^ ^d+2-2n-ne _ rj.^^ j^p2 
equation is then one-loop ultraviolet renormalizable for 
d < 4 + 39. Thus, the one-loop ultraviolet renormaliz- 
ability can be extended to dimensions four and greater 
depending on how the noise scales in the ultraviolet re- 
gion. We will not pursue these issues further in this paper 
and henceforth restrict our attention to white noise. 

There is a formal connection with X{(j)'^)[d,+i]+i QFT 
that is made apparent by extracting a factor of v^k^ from 
equation (p^. If we do so, we obtain 

VM = 




k^ + ^fA+ 0{A^). 



(17) 



This is recognizable as the effective potential for A^^ 
QFT in d -|- 1 spatial dimensions, that is, Euclidean 
[d + 1] + 1 spacetime dimensions |l^Jl^. To make the 
connection, interpret XFq/v'^ as the mass term of the 
QFT, \^ /v"^ as Aqft (the coupling constant of the QFT), 
and V as the field 0qft- 

After renormalization the tadpole is to be set to zero, 
but we still need a bare tadpole to act as a counterterm 
during the intermediate stages of the regularization. This 
is completely analogous to the situation in massless 0^ 
QFT, where the renormalized mass of the quantum field 
is fine tuned to zero by hand [|l|,0. 

We now evaluate the effective potential in different spa- 
tial dimensions. 

Case d= \: We are interested in the integral 



I{a) = 



dfc2 



(18) 



This integral is divergent and needs to be regularized. 
The most direct way to do so is by the "differentiate and 
integrate" trick which leads to [Oj 
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T{a) = na — —a 



(19) 



Here k is an infinite constant of integration. We absorb 
K into the bare action, where it renormalizes Fq [ p^jis) 



V\v:d- 



[{Fo 



1 A3/2 
T-A^r- 



-0{A^). 



{Fo + iA«2) 



2\3/2 



F 



3/2 



(20) 



These are all renormalized parameters at 0{A). The 
zero-loop contribution will always dominate at large 
fields (rapid fiow). Near w = 0, it is the one-loop contri- 
bution that is dominant. We take the physical value of 
the renormalized tadpole, Fq = 0, and encounter some- 
thing very interesting - the system undergoes dynamical 
symmetry breaking (DSB) in a manner qualitatively sim- 
ilar to the Coleman- Weinberg mechanism of QFT. The 
effective potential then simplifies to 



|wp + 0(y^2). (21) 



The potential is not analytic at zero field (a phenomenon 
well known from massless QFTs) . For large v the 

classical potential dominates, whereas for small v one- 
loop effects dominate. Thus the symmetric vacuum (v = 
0) is unstable. While it is easy to see that the symmetry is 
dynamically broken, the presence of the unknown 0{A^) 
terms make it difficult to make a quantitative estimate 
for Wniin- We can nevertheless qualitatively estimate the 
shift in the expectation value of the velocity field 
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A 



27r 21/2 j,2 



0{A^). 



(22) 



This DSB is particularly intriguing in that it suggests 
the possibility of a noise driven pump. For example, in 
thin pipes (where the flow is essentially one-dimensional) 
and provided the physical situation justifies the use of 
Burgers equation, this result indicates the presence of a 
bimodal instability leading to the onset of a fluid flow 
with velocity depending on noise amplitude. 

Case d = 2: Because of the analogy between the one- 
loop effective potential for the KPZ equation (for white 
noise) and that for the Xcf)"^ QFT, we can write down 
the renormalized one-loop effective potential by inspec- 
tion il4|,|Tl 



V[v;d 




\{^l)v^f-[F,{^,)f] 
[Fo(/i) + kmv^f 



X lOE 



[Fo(Ai)]2log 



i^o(/i) 



+0{A^). 



(23) 



Here /i is the renormalization scale in the sense it is used 
inQFT|l|,0. 

We now tune the renormalized tadpole to its physical 
value of zero. After performing a flnite renormalization 
to simplify the expression, we obtain 



\2 



1 



(27r)2 81/ 



V log ^ 



0{A^). 
(24) 



The one-loop contribution will always dominate at large 
fields (rapid flow). Near v = Q the logarithmic singular- 
ity is rendered finite by the polynomial prefactor. That 
the KPZ potential has a non-trivial minimum exhibit- 
ing DSB is exactly the analog of the Coleman- Weinberg 
mechanism la. We can estimate the location of the 
minimum by differentiating the effective potential. We 
obtain 



±/j, exp 



2 A2 ^ 4 + 



(25) 



Note that Umin ^ as ^ 0, to recover the tree level 
minimum Wmin = 0. The present discussion is relevant 
to either (1) surface evolution on a two dimensional sub- 
strate, or (2) thin superffuid films, since superfiuids are 
automatically vorticity free, justifying the application of 
the KPZ equation. 

We can also calculate the one-loop beta function, for 
as the bare effective potential does not depend on the 
renormalization scale, one has that /idV/d/i = 0. We 
get 



„ dof d , ^ A^ 2 
^^ = ^d^'=4^^' + ^^-^) 



(26) 



Case d = 3: We are interested in the integral 



2 \ dcf 



dfc^fc^ 



The "differentiate and integrate" trick yields 

'I{a?) = nia? + ^20^ H . 
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(27) 



(28) 



Here ki and K2 are two infinite constants of integration. 
Since the sign in front of the term is positive there is 
no possibility of DSB in d = 3. We absorb ki and K2 into 
the bare action, where they renormalize both and A 



V[v-d = ?,]=^^[{Fo + \\v^f-F^] 



(30^2 

-0{A^] 



-A- 



A5/2 



,5/2 



-F, 



5/2 



(29) 



These are all renormalized parameters at order O(^). 
For large fields the one- loop correction is dominant, 
whereas for small fields the tree level result dominates. 
The effective potential is positive and monotonically in- 
creasing. 

If we set the renormalized tadpole to its physical value 
of zero, we have 



\2 

V[i;;d = 3] = —V 



1 



A- 



A5 



(307r2) 25/2;/4 



-\vf + 0{A^). (30) 



At zero- loops the vacuum is symmetric at Wmin = 0. One- 
loop physics does not change this. There is an important 
sign change in the one-loop contribution when comparing 
d = 3 with d = 1. The DSB that is so interesting in 
d = 1 and d = 2, is now absent in d = 3. Notice that the 
effective potential is non-analytic at Wmin = 0. This case 
is of most interest in hydrodynamics and we have shown 
that (at one-loop) the ground state (wmin = 0) is stable 
under Gaussian noise. 

Discussion: When taken as a model for turbulence, it 
is natural to study the ultraviolet properties of the KPZ 
equation because in a turbulent fluid, energy cascades 
down from larger to smaller length scales and one might 
expect scale-invariant solutions in the limit as A: ^ -l-oo. 
Indeed, Forster, Nelson, and Stephen pointed out long 
ago the importance of studying the short distance and 
short time (i.e., ultraviolet) correlations in the fluid ve- 
locity in fully developed turbulence jj]. More recently, 
Yakhot [0 undertook an analysis of the ultraviolet prop- 
erties of the randomly forced Navier-Stokes (NS) equa- 
tion and found that ultraviolet renormalizability requires 
modifying the fluid dynamical equation by a term not 
originally 'present in the NS equation, though this new 
term is much more complicated than the tadpole needed 
to successfully renormalize the KPZ equation. 

The ultraviolet properties of the KPZ equation ex- 
posed in this paper reveal a rich structure and comple- 
ments studies concerned with the long distance (infrared) 
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features of this equation. This distinction shows up in a 
striking way in both, the need of a bare tadpole for ul- 
traviolet renormalizability and the presence of dimension 
dependent dynamical symmetry breaking for the configu- 
rations studied here. Our focus in this paper has been on 
the concept of the effective potential and we have shown 
how this concept, so useful in QFT, can be extended to 
stochastic equations containing gradient forcing terms, 
such as the KPZ equation, by employing inhomogeneous 
and static fields. 

The evidence for the dynamical symmetry breaking 
found here must be interpreted with due caution. We 
must emphasize, for the reasons explained earlier, that 
the effective potential is calculated for a restricted class 
of field configurations and the symmetry breaking is with 
respect to this class of constant fluid velocity configura- 
tions. Thus, while for d = 3 the one-loop effective po- 
tential is minimized for a zero-velocity background flow 
^min = 0; the potential in both d — 1 and d = 2 is 
minimized at a non-zero value of the fluid flow, namely 
( p^ ) and (|2^), respectively. The symmetry breaking is 
thus manifested as a preference for the system to "jump" 
from zero to a non-zero but constant flow. The jump is 
triggered, of course, by the noise source. In the surface 
growth interpretation, the minimum of the potential for 
d = 3 corresponds to a surface with zero slope, while the 
dynamical symmetry breaking that occurs in d = 2 and 
d — 1 corresponds to a jump in the surface slope to a 
constant non-zero value. This should be contrasted to 
recent exact and numerical calculations of the PDF for 
the KPZ equation in d = 1 dimension. By making use 
of a lattice version of the KPZ equation in one dimen- 
sion, Derrida and Lebowitz have performed an exact cal- 
culation of the probability distribution of the deviation 
of the average current Q . The results of these authors 
agree with a numerical calculation of the KPZ PDF in 
d = 1 by Yakhot and Chekhlov In these papers, 

the PDF, which is essentially the exponential of the ef- 
fective potential, gives information about the probability 
to find a velocity difference Aur = v{x + r) — v{x) at two 
widely separated points in the fluid (i.e., as r ^ -l-oo) 
and for asymptotically large times. Their PDFs have a 
global maximum at Au^. — 0, corresponding to no ve- 
locity difference, i.e., a spatially constant velocity (or a 
constant slope surface). Thus, the most probable config- 
uration corresponds to a spatially constant velocity. Of 
course, these PDFs take into account all dynamically al- 
lowed configurations that evolve via the KPZ equation; 
our PDF, by contrast, samples only the class of constant 
velocity configurations. Nevertheless, and within these 
limitations, we have been able to show that for the lin- 
ear theory (A = 0) the ground state has zero velocity (or 
zero slope) and that once we turn on the nonlinear inter- 
action, and restrict the fluid to have a constant velocity, 
there is a dynamically preferred velocity as determined 
by one-loop physics in d = 1 and d = 2 dimensions. 

In this paper we have also demonstrated that the KPZ 
equation with white noise is at least one-loop renormaliz- 



able for d = 1,2,3 dimensions but is non-renormalizable 
(at one-loop) for d > 4. In making statements about 
renormalizability, care must be taken in distinguishing 
ultraviolet and infrared properties. Thus, from the work 
of Sun and Plischke ||l6| and Frey and Tauber ||l^, it 
is known that the KPZ equation is, at least, two-loop 
renormalizable for d = 1,2,3. Although these authors 
were interested in the infrared properties of the KPZ 
equation, they had to deal with both ultraviolet and in- 
frared divergences in their work. By making use of the 
mapping of the KPZ equation to directed polymers with 
quenched noise, the one-loop renormalizability has been 
demonstrated in 2 < d < 4 . 

Though the calculations presented here are limited to 
one-loop order in the noise amplitude, it must be em- 
phasized that in many cases one-loop physics is enough 
to capture essential features of stochastic partial differ- 
ential equations (SPDEs) |0 and QFTs. Furthermore, 
while calculations of the one-loop effective potential are 
straightforward, one should not forget that it is neces- 
sary to develop a physical interpretation of this effective 
potential to ensure that it is as relevant for SPDEs, as 
it is for QFTs |10|. In this paper we have presented an 
analogy between the statics of the KPZ equation and the 
static behavior of Xcj)'^ QFT -the vacuum structure of the 
(f>* QFT carries over into the basic field configurations of 
the KPZ equation. In 1 and 2 space dimensions we have 
exhibited the occurrence of DSB, which does not take 
place in 3 space dimensions. 
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